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ABSTRACT 

We propose a new method to measure the metalhcity of Type la supernova progenitors using Mn 
and Cr hues in the X-ray spectra of young supernova remnants. We show that the Mn to Cr mass 
ratio in Type la supernova ejecta is tightly correlated with the initial metallicity of the progenitor, 
as determined by the neutron excess of the white dwarf material before thermonuclear runaway. We 
use this correlation, together with the fl ux of the Cr and Mn K g X-ray lines in the Tycho supernova 
remnant recently detected by Suzaku (jTamagawa et al.l l2008l ) to derive a metallicity of log{Z) ~ 



1.32„Q 33 for the progenitor of this supernova, which corresponds to log{Z/Z(i 



0.60 



+0.31 
-0.60 



according 



to the latest determination of the solar metallicity bv lAsplund et all (j2005f ). The uncertainty in the 
measurement is large, but metallicities much smaller than the solar value can be confidently discarded. 
We discuss the implications of this result for future research on Type la supernova progenitors. 
Subject headings: stars:binaries:close — supernova remnants — supernovaergeneral — X-rays:ISM 
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1. INTRODUCTION 

Despite decades of continuing effort, the nature of the 
progenitor systems of Type la supernovae (SNe) remains 
unknown. The most widely accepted theoretical scenar- 
ios involve the thermonuclear explosion of a C+0 white 
dwarf (WD) destabilized by accretion of material from 
a binary companion, cither another WD (double degen- 
erate systems, DDs) or a normal star (single degener- 
ate systems, SDs). Observations of supernova rates in 
distant galaxies sugge st at least two dif ferent wajs to 
produce Type la SNe (jScannapieco fc BTldst cn 20051 ): a 
'prompt' channel associated with young stellar popula- 
tions and a 'delayed' channel associated with old stellar 
populations. At present, it is not clear what relationship, 
if any, these channels have with the theoretical scenar- 
ios, or even if the propose d scenarios c an yield Type la 
SNe at the observed rate (|Mao j l2008l ) . Most attempts 
to constrain the fundamental properties of the progeni- 
tors have been unsuccessful, includi ng direct identifica- 
tion in optical pre-explosion images aoz fc Mannuccil 
|2008|) , searches for H stripped from the binary companion 
in SN spectra ()Leonardll2007f ). and searc hes for prompt 
emiss i on from circumstel lar interaction (jPanagia et al.l 
120061 iHughes et aIll2007D . Very recently, the possible 
identification of the progenito r in a pre-exp losion X-ray 
image of SN 2007on by Voss fc Nelemansj (120081) has; 
been questioned by post-explosion images |Roelofs et "all 
|2008() . Campaigns to look for the surviving companion 
star to the exploded WD in the nearby Tycho supernova 
remnant (SNR), known to be of Type la origin, have 
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also p roduced controversial results (jRuiz-Lapuente et al.l 
l2004tllhara et al.ll2007t ). This lack of clues to the identity 
of the progenitors in the observations of Type la SNe is 
sometimes referred to as 'stellar amnesia'. 

One of the most important constraints on the age and 
nature of the progenitor systems is their metallicity Z. 
During stellar evolution, the C, N, and O that act as 
catalysts for the CNO cycle pile up into ^"^N, which is 
converted to ^^Ne in the hydrostatic He burning phase 
through the reactions '^'^'N{a,j)^^F{f3+ ,Ve)'^^0{a,j)'^'^'Ne. 
The Z?"*" decay of ^^F increases the neutron excess of the 
WD material, defined as 77 = 1 — 2{Za/A) (where Za is 
the atomic number and A is the mass number), resulting 
in a linear scaling of tj with metallicity: rj = 0.101 x 
z (jr immes et all 120031 ) ■ This can have an observable 
impact on Type la supernova spectra (jLentz et al.ll2000( ). 
but the predicted effects are hard to detect in practice (in 
one case, SN 2005bl, a progenitor of subsolar m etallicity 
has been suggested bv lTaubenberger et al.|[2"008f ). In this 
Letter, we describe a new direct method to measure the 
metallicity of Type la SN progenitors using Mn and Cr 
lines from shocked ejecta in the X-ray spectra of SNRs. 

2. THE Mmn/Mcr ratio AS A TRACER OF 
PROGENITOR METALLICITY 

2.1. Mn/Cr vs. Z 

Nuclear burning in Type la SNe happens in different 
regimes depending on the peak temperature reached by 
the WD material. From higher to lower temperatures, 
the regimes are: nuclear statistical equilibrium (NSE), 
incomplete Si burnin g, incomplete O burning and incom- 
plete C-Ne burning (|Wooslevlll986[ ). During the explo- 
sion itself, electron captures are too slow to change the 
original value of 77 except in the innermost ~ 0.2 Mq of 
ejecta, where nuclear burning h appens at higher dens i- 
ties, leading to neutron-rich NSE (jBrachwitz et al.l2000( l. 
Outside this region, an imbalance between free neutrons 
and protons docs not affect the production of the most 
abundant species, with t he possible exceptio n of ^^Ni at 
supersolar metallicities (jTimmes et al. 1 120031 ). However, 
some trace nuclei with unequal numbers of protons and 
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Fig. 1 . — Mn to Cr mass ratio as a f unctio n of progenitor metallic- 
ity for Type la IIBadenes et al.ll2003l . [20051. this work) and core col- 
lapse llWooslev fc We^^iTIT995l) SN models, represented with red 
circles and blue triangles, respectively. The red and blue shaded 
regions encompass all the models in each class. The dashed line 
represents a power-law fit to the Type la models. The orange and 
green plots illustrate the impact of weak interactions during the 
simmering phase for very subluminous Typo la SNe (Amis = 1-8) 
in WDs with main seq uence progenitor masses of 3 Mq and 6 Mq, 
respectively (sec ^ 12.21 1. The horizontal solid and dotted lines rep- 
resent the estimated Mm^/Mcr ratio in the Tycho SNR from the 
Suzaku observations (0.74 ± 0.47, shaded in gray), which corre- 



sponds to log(Z) 



-1.32 



-1-0.67 



(see §1211. 



neutrons are very efficient at storing the neutron excess. 
The most abundant among these nuclei is ^^Mn, which is 
produced during incomplete Si burning as ^^Co. When 
normalized to another product of incomplete Si burning 
whose nucleosynthesis is insensitive to the neutron ex- 
cess (Cr is the ideal choice), the yield of Mn becomes an 
excellent tracer of the progenitor metallicity. 

This diagnostic can be quantified with Type la SN 
explosion mod e ls. We hav e used the model grid from 
I Badenes et al.l ()2003L l2005| ). complemente d with other 
model s calculated with the code described in lBravo et al.l 
(Il996l). Physica l input s and initial conditions are as 
in iBadenes et al.l ()2003[ ). except in the cases noted be- 
low. We consider examples of several explosion mecha- 
nisms: deflagrations, delayed detonations, and pulsating 
delayed detonations. For the present work, we have re- 
calculated the nucleosynthesis of four delayed detonation 
models (DDTa, DDTc, DDTe, and DDTf) at different 
metallicitics by altering the value of X(^^Ne) in the prc- 
explosion WD from the canonical 0.01 (corresponding to 
Z = 9.0 X 10"'^). We have also calculated one deflagra- 
tion model starting from a central density twice as large 
as the other models, and one delayed detonation model 
with the C to O mass ratio Mc /Mq = 1/3 in the WD 
instead of the usual 1. In spite of these fundamental dif- 
ferences in explosion mechanisms and initial conditions, 
we find a very tight correlation between the Mn to Cr 
mass ratio outside the neutron-rich NSE region in the 
SN ejecta, MMnl^Cr, and the progenitor metallicity Z 
(see Figure [T]). A power-law fit yields the following rela- 
tion (with a correlation coefficient = 0.9975): 



MMn/Mcr = 5.3 X Z 



0.65 



(1) 



Removal of the inner 0.2 M0 of neutron-rich NSE ma- 
terial from this relation is justified because (1) there 
is observational evidence that this material does not 




Fig. 2. — Upper limits to M^onv during the simmering phase 
(left axis, black plots) and its overlap with the Si rich region of 
ejecta (right axis, colored plots) as a function of Z. The black 
plots represent quadratic interpolations for different Mms'- 3Mq 
(dashed), 5Mq (dotte d), and GMp) (solid), wi th the symbols indi- 
cating the values from lDommguez et al.l 1)20011 ) (for Mms = 6Mq 
at Z = 0.02, we have taken the average between Mms = 7M0 
and Mms = SMq). The colored plots represent the overlap (in 
%) between M^onv and the Si-rich region of ejecta, for very sub- 
luminous (Amis = 1.8, red), mildly subluminous (Amis = 1.6, 
orange), and normal, but faint (Amis = 1.4, blue) Type la SNe. 



mix outwar ds during the explosion (jCcrar dv et al."2007|; 
lMazzali"era l. 2007) or its aftermath (Fcscn et a l. 2007i) ; 
and (2) the dynamical ages of most ejecta-dominated 
Type la SNRs are too small for the reverse shock to have 
reached so deep into the ejecta. For illustration, we have 
plotted in Figure |l ] the core collapse SN models from 
IWooslev fc Weaveil ()1995f ). which also show some corre- 
lation between MMn/Mcr and Z, albeit with a much 
larger spread and a much shallower slope. 

2.2. The impact of C simmering 

In order for Eq. [T]to hold, the value of r] set by the pro- 
genitor metallicity must remain unchanged between the 
formation of the WD and the SN explosion. This should 
be true for DD progenitors if the final merger and run- 
away happen on dynamical timescales. In slowly accret- 
ing WDs, 77 can be modified through electron captures 
during the so-called 'simmering' phase of non-explosive C 
burning that tak es place in the 1000 yr prior to the ex- 
plosion (jPiro fc Bildsten 2008.). The impact of this addi- 
tional neutronization on the Mmu/ Mcr ratio will depend 
on the extent of the convective region over which the neu- 
tronized material is mixed {Mconv), and its overlap with 
the portion of the WD that will undergo incomplete Si 
burning where Mn and Cr are sy nthesized. 

Some previ ous st udies llKuhlen et all l2006l : 
iPiro fc BildstCTi I2OO8I : IChamulak et all I2008D have 
found large values for Mconv by assuming an homo- 
geneous chemical composition for the WD, effectively 
applying the Schwarzschild criterion for convection. We 
will adopt the hypothesis that the convective region 
is limited by the Ledoux criterion to the C-depleted 
core of the W D created during hyd rostatic He-shell 
burning, Mcore ijHoflich fc Steinll2002D . In the absence 
of a self-consistent picture for con vection inside WDs , 
this must remain an open issue (see[P iro fc Chang|[2008l . 
for a discussion), but we believe that our hypothesis is 
reasonable given the existence of other processes that 
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Fig. 3. — X-ray spectrum of the Tycho SNR observed by Suzaku 
in the vicinity of the Fe Ka line. The Ka Hnes of Cr and Mn 
are detected at > lOcr and 7(t confidence levels, with fluxes of 
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tively (|Tamagawa et al.ll2008h . 

can reduce the extent of Mconv and limit the mixing 
of neutroni zed material. These in clude the presence of 
Urea shells (IStein fc Wheelei1l2006f ) and the long thermal 
diffusion time scale for buoyant bub bles larger than 
- 100 m jG arcia-Senz fc WooslevlfT995l ) . which prevents 
them from mixing completely with their surroundings 

during convectiom 

According to iMazzali et al.l ()2007[ ). Si-rich ejecta in 
Type la SNe lie between Lagrangian mass coordinates 
1.05 Mq and (1.55 - 0.69 x Atoi5)Mq, where Amis 
is the light curve width parameter (defined as the de- 
cline in blue magnitude 15 days after maximum), which 
acts as a proxy for SN brightness. In this context, the 
extent of the overlap between Mconv and the Si-rich 
ejecta depends on the main sequence mass M]\js and 
initial meta Uicity of the WD prog enitor (which deter- 
mine Mcore, iDomi'nguez et al.ll2001l ) and the SN bright- 
ness (see Figure [2]) . This overlap is only significant for 
subluminous (Atoi5 > 1.6) Type la SNe originated by 
progenito rs with either large Mms or low or both. 
Oiamulakjrn] dlOOS) find an upper limit for the in- 
crease of rj during the simmering phase of A77 = 0.0015, 
which is comparable to the value of rj in solar ma- 
terial. The impact of simmering on the MMn/Mcr 
ratio can then be estimated by mixing material with 
MMn/M cr = 0-3 (appropriat e for the value of Zq de- 
rived by lAsplund et al.l |2005[ ) into the incomplete Si 
burning region, in a proportion equivalent to the extent 
of the overlap shown in Figure [21 The green and orange 
plots in Figure [T] are two examples of such 'simmering- 
modified'' models for very subluminous (Atoi5 = 1.8) 
Type la SNe, illustrating our conclusion that C simmer- 
ing will only modify the Mmu/Mct ratio in the SN ejecta 
for very subluminous SNe, and then only in cases where 
Mms is large, or Z is low, or both. 

3. MEASURING THE Mmn/Mcr RATIO IN SNRS 

The work in this Letter is motivated by the recent 
Suzaku detection of Mn and C r in the X-ray sp e ctrurn 
of the Tycho SN reported by iTamagawa et al.l (|2008f ) 
(see Figure |3|). This observational result opens the 
possibility of studying the M]\jn/Mcr ratio in Type 
la SN ejecta, which cannot be done using optical SN 
spectra due to the 2.7 yr half- life of ^^Fc in the de- 
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Fig. 4. — Interpolated Mn to Cr specific Ka emissivity ratio as 
a function of rigt and kT . The dotted box encompasses the val- 
ues of riet and kT found in the Si-rich ejecta of six young Type 
la SNRs bv lBadenes et al.l ll2007h : 9.49 < login^t) < 11.94 cm'^s; 
1.0 < kT < 10.0 keV. The striped area corresponds to the re- 
gion of parameter space appropriate for the Tycho SNR, 10.23 < 
log{net) < 10.99 cm-^s; 1.0 < kT < 10.0 keV. 

cay chain ^^Co-^^^Fe^'''^Mn. Since Mn and Cr are 
synthesized together in the explosion and have very 
similar electronic structures, it is possible to estimate 
their mass ratio from the line flux ratio: Mmn/Mcr = 
1.057 X {FmuI Fcr)/ {EmuI Ecr), where 1.057 is the ra- 
tio of atomic masses, Fmu/ Fcr is the line flux ratio, and 
EMn/Ecr is the ratio of specific emissivities per ion. 

Public data bases for X-ray astronomy do not usu- 
ally include lines from trace elements like Mn and Cr, 
but the value of E]\in/Ecr can be estimated by inter- 
polation along the atomi c number sequence from ele- 
ments with available data (jHwang et al.ll2000l). W e have 
used the ATOMDB data base (jSmith et al.ll200lD to re- 
trieve Ka line emissivities for Si, S, Ar, Ca, Fe, and 
Ni (Za = 14, 16, 18, 20, 26, 28) as a function of ioniza- 
tion timescale n^t and electron temperature kT, the two 
variables that control the line emission of a plasma in 
noncquilibrium ionization. Then we performed a spline 
interpolation to obtain Emu/ Ecr, which we plot in Fig- 
urelH together with the region of the {net, kT) parameter 
space that is populated by young Type la SNRs (see the 
Appendix in iBadenes et al.ll2007t ). The uncertainty in 
the value of Emu/ Ecr comes from the variation within 
this region and the error introduced by the interpolation 
itself, which we have estimated by comparing the true 
and interpolated values for the Eat/Es ratio. For the 
larger region of the {uet, kT) plane plotted in Figure IH 
we find Emu/Ect = 0.69±0.32 ( Emu/Ect = 0.66±0.26 
for the smaller region appropriate for Tycho). 

4. DISCUSSION AND CONCLUSIONS: THE METALLICITY 
OF TYCHO'S PROGENITOR 

From the Suzaku observation, FMn/Fcr = 0.46 ±0.21, 
which gives Mun/Mcr = 0.74 ± 0.47. This mass ra- 
tio translates into a metaUicity of Z = 0.048^Q g3g for 

the progenitor of the Tycho SNR {log{Z) = -1.32l° j^). 
The large error bar in this result is dominated by the sta- 
tistical uncertainties in the X-ray fluxes from the faint 
Mn and Cr lines, and it will be reduced in an upcom- 
ing, approved deeper Suzaku observation. At present, 
we find a strong indication for a supersolar metallic- 
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ity {log{Z/ZQ) = 0.45 
iGrevesse fc Saurol ([l99l); log^Z/Zo) 



Q with the solar value from 



0.60 



+0.31 
-0.60 



with 



the newer value from lAsplund et"al] (|2005[ )). Our results 
are compatible with a solar metallicity, but subsolar val- 
ues can be discarded with confidence. Given this mea- 
surement and the fact that Tycho's SN was probably 
eithe r normal or slightly overluminous (iBadenes et al.l 
I2006t iRuiz-Lapuentelbool . we do not have to concern 
ourselves with the impact of neutronization during C 
simmering in this particular case. The Tycho SNR is 
59 pc above the Galactic plane at a Galactoc entric radius 
of 9.4 kpc (assuming a distance of 2.4 kpc, ISmith et al.l 
Il991t) . A supersolar metallicity is higher than average 
for this locatio n, but well within the sp read of measured 
[Fe/H] values (|Nordstr5m et all [2001 ). This combina- 
tion of Galactrocentric radius and metallicity suggests a 
young progenitor age (a few Gyr or less), which would 
make Tycho a candidate for the 'prompt' channel of Type 
la SNe, but the scatter in the age-metallicity relations 
and the uncertainties in our measurement are too large 
to be more specific on this point. 

In this Letter^ we have proposed a new method to mea- 
sure the metallicity of Type la SN progenitors using Mn 
and Cr lines in the X-ray spectra of their SNRs. We 
have applied it to the Tycho SNR and obtained a strong 
indication that its progenitor had a solar or supersolar 
metallicity. The main strength of our method is its sim- 
plicity: it is based on well-known nuclear physics and 
observable parameters that are easy to measure. Detec- 
tion of Mn and Cr lines in other young Type la SNRs 
in the Galaxy should be possible with Suzaku, Chan- 
dra, and XMM-Newton observations. A few additional 
SNRs in the Magellanic Clouds might be reached with 
deep exposures, bringing the total of potential targets 
to perhaps a dozen. This number would increase sig- 



nificantly with the inclusion of SNRs in nearby galax- 
ies like M31 and M33, which may be accessible to next 
generation X-ray observatories like Constellation- X and 
XEUS. Even with a small sample of objects, we should 
be able to use this method to verify and complement the 
indirect metallicity studies of extragalac tic Type la SNe 
(jCallagher et al.ll2005t[Prieto et aLl^OOl ) and test claims 
about the metallicity d ependence of the Type la SN rate 
(|Kobavashi et al.lll998l ). Measurements or upper hmits 
below a certain threshold {Almn/Mcr ^ 0.1) would also 
provide interesting constraints on the extent of the con- 
vective region in accreting WDs. We conclude by not- 
ing that Mn and Cr lines have alrea dy been detected b y 
ASCA in the Galactic SNR W49B (jHwang et al.|[2000l ). 
The measured line flux ratio also suggests a solar or su- 
persolar progenitor metallicity, but b oth the age and the 
SN type of W49B a re controversial (jHwang et al.ll2000l : 
IBadenes et al.l [20071) . We defer a detailed discussion of 
this object and the application of our method to other 
Galactic SNRs to a forthcoming publication. 
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